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ABSTRACT 

We present ultraviolet, optical, near- infrared, Spitzer mid-infrared, and radio 
images of 14 radio-selected objects in M 33. These objects are thought to repre- 
sent the youngest phase of star cluster formation. We have detected the majority 
of cluster candidates in M 33 at all wavelengths. From the near-IR images, we 
derived ages 2-10 Myr, i^ 5 -band extinctions (A Ks ) of 0-1 mag, and stellar masses 
of 10 3 -10 4 M Q . We have generated spectral energy distributions (SEDs) of each 
cluster from 0.1 /im to 160 /im. From these SEDs, we have modeled the dust 
emission around these star clusters to determine the dust masses (1-10 3 M Q ) and 
temperatures (40-90 K) of the clusters' local interstellar medium. Extinctions 
derived from the JHK$, Ha, and UV images are similar to within a factor of 
2 or 3. These results suggest that eleven of the fourteen radio-selected objects 
are optically-visible young star clusters with a surrounding H II region, that two 
are background objects, possibly AGN, and that one is a Wolf-Rayet star with a 
surrounding H II region. 

Subject headings: galaxies: individual (M 33) — galaxies: star clusters — ISM: 
dust, extinction — ISM: structure — optical: ISM — infrared: ISM 

1. INTRODUCTION 

A galaxy's morphological evolution and the properties of its interstellar medium (ISM) 
are driven by massive stars in young (0-10 Myr) star clusters. At these ages, star clusters 
produce metals from supernovae explosions, generate ionizing photons that determine the 
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phase of the ISM, and provide mechanical energy to rearrange the ISM. These effects originate 
in the local ISM but eventually contribute to galaxy-wide ISM changes (see Calzetti et 
al. 2004). Thus, the small-scale effects of young star clusters on their ISM have evolutionary 
consequences on the host galaxy in general. 

When a young star cluster first forms, the OB stars in that cluster have not had time to 
remove their natal molecular cloud material, either by photoionization or by the mechanical 
force from the stellar winds. Because this molecular material exhibits optical depths greater 
than 1 at visible wavelengths, the cluster should be detectable only at infrared (IR) and radio 
wavelengths. Johnson et al. (2001) found 14 thermally-emitting objects in M 33 using high 
spatial resolution radio continuum maps at 6 cm and 20 cm. These selected candidates have 
a spectral index (a^^J °^ ^ or S rea t er > indicating optically thick free-free emission from 
a thermal source. The Lyman continuum photon rates determined from radio continuum 
maps imply the presence of 2-4 07 equivalent stars and hence an age of 7 Myr or less. 

To learn more about these radio-selected objects, we undertook a multiwavelength sur- 
vey from the radio to the far- ultraviolet (far-UV). With the combination of Spitzer InfraRed 
Array Camera (IRAC; Fazio et al. 2004) and Multiband Imaging Photometer for Spitzer 
(MIPS; Rieke et al. 2004) images, we identify the radio-selected objects from Johnson et 
al. (2001) and produce spectral energy distributions (SEDs) to determine their dust proper- 
ties (temperature and mass) as well as morphologies of the polycyclic aromatic hydrocarbon 
(PAH) /dust emission. With a near-IR color-color plot and absolute magnitudes, we deter- 
mine their stellar ages, masses, and extinctions. Finally, archival Ha: and ultraviolet (UV; 
GALEX) images were used to detect objects at shorter wavelengths for which we measured 
UV and Ha extinctions based on stellar synthesis models. The observations and data reduc- 
tion are presented in §2. Basic results of our analysis are given in §3, and discussion of our 
findings is given in §4. In §5 we summarize our findings. 



2. OBSERVATIONS AND REDUCTIONS 

2.1. Spitzer IRAC and MIPS Images 

Spitzer MIPS and IRAC images of M 33 were produced under the Guaranteed Time 
Observer (GTO) program "M 33 Mapping and Spectroscopy" (PID 5) by Spitzer Science 
Working Group member R. D. Gehrz. Hinz et al. (2004) describe the MIPS observations 
used in this article. Observations of M 33 were made using all four bands of IRAC on 
2004 January 9. The M33 IRAC mapping sequence consisted of 438 frames per channel, 
including a 3 point 1/2 pixel dither for each map position. The integration time was 12 s 
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per frame. The raw Spitzer IRAC data (single epoch) were processed and flux calibrated 
with version 10.0.1 of the Spitzer Science Center (SSC) pipeline. Post-basic calibrated data 
(Post-BCD) processing was carried out using the 2004 June 6 Linux version of the SSC 
Mosaicking and Point source Extraction (MOPEX) software (Makovoz & Khan 2005). Four 
steps of MOPEX were implemented: cosmetic fix, background matching, outlier detection, 
and mosaicking. The cosmetic fix module reduced column pull-down and mux-bleed (trails 
that form rightward of bright objects) artifacts associated with 3.6 /im and 4.5 /im data. 
Background matching was performed by minimizing the pixel differences in overlapping areas 
with respect to a constant offset computed by the program. Cosmic rays and other outliers 
were detected and eliminated with the outlier detection module. In the final step, the images 
were reinterpolated to a pixel scale of approximately 0'/86 pixel -1 . The final science images 
are in units of MJy sr _1 . 



2.2. SQIID JHK S Images 

The near-IR JHK S data were obtained 2002 October 16-23 on the 2.1 m telescope at 
Kitt Peak National Observatory using the NOAO SQIID near-IR instrument. SQIID obtains 
simultaneous images in JHK$ at a pixel scale of 0'.'7 pixel -1 over a 304 " x 317 " area. One 
night was photometric while another night was sporadically photometric. We observed each 
radio-selected object position with a standard 5 s integration time, 24 coadds, and 3 sets of 
5 dithered images (total integration times of 75 s). Throughout these nights, we observed 
3-6 UKIRT faint standard stars (Hawarden et al. 2001) for flux calibration. 

All data reductions were performed using the Image Reduction and Analysis Facility 1 
(IRAF) software. Averaged dark images were subtracted from individual images, followed by 
flatfielding with the averaged flatfield image for the appropriate filter. We created sky images 
from each set of dithered science images and then subtracted the scaled sky background 
images from each science image in their respective sets. The science images were then 
registered and combined. Finally, a world coordinate system was generated for each JHK S 
image from field stars' celestial coordinates in the 2MASS point source catalog. 

Flux calibration was performed with the standard stars. Zero points were determined 
using a 7" aperture; we found that color terms were negligible. Between the flux calibration of 
the photometric and near-photometric nights, the zero point differences were 0.01 magnitudes 



1 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the As- 
sociation of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 
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for J and H, and 0.08 magnitudes for K s . We compared our derived magnitudes of field 
stars to those found in 2MASS, appropriately transformed to the UKIRT system (Cutri et 
al. 2003). These magnitudes were found to be in close agreement (~0.05 mag for all filters). 

2.3. Archival Imagery 

We searched multiwavelength data archives for images of these radio objects. We found 
four other wavelength bands that had suitable data. VLA 6 cm continuum maps of the 
fourteen radio-selected objects were generated from the data taken for Duric et al. (1993). 
These radio maps were used to discern which object in the Spitzer images is the actual 
compact radio object, and these radio maps were also used by Johnson et al. (2001). The 
SQIID near-IR observations did not cover all M 33 radio-selected objects. Therefore, we 
used the M 33 2MASS mosaic produced by Jarrett et al. (2003). This image allowed us to 
derive JHK S magnitudes and colors for Johnson et al. (2001) objects 1-4 and 12. R-band 
and Ha images of M 33 were taken from Massey et al. (2002). Continuum emission in the 
Ha image was removed by subtracting a scaled R-band image. We determined the scaling 
factor by comparing the aperture fluxes of 20 stars in both images. To determine what the 
corresponding flux is for one count (data number or DN) in the continuum-subtracted Ha 
image, we used Ha fluxes of three round, clearly separated H II regions given in Hodge et 
al. (2002, ; BCLMP 2 71a, 95 Ab, and 633a+b). Hodge et al. (2002) give Ha fluxes as well as 
the areas for the H II regions. We choose isolated, spherical H II regions in order to both avoid 
contamination from other sources and to make use of simple circular aperture photometry. 
The conversion factor derived is 7.64xl0~ 19 ergs s _1 cm~ 2 DN -1 . Far-UV (0.1528 /im) and 
near-ultraviolet (near-UV; 0.2271 /im) archival images of M 33 were downloaded from the 
GALEX data archive to be used for aperture photometry. Flux calibration of the GALEX 
data from DN s _1 to Jy is given in the GALEX data handbook 3 . 

3. RESULTS 

3.1. Finding the Radio-Selected Objects 

We located the 14 objects from Johnson et al. (2001) in the far-UV (0.1528 /im), near- 
UV (0.2271 /im), Ha, JHK S) and Spitzer images. Figures 1 through 14 show 90" sqaure 



2 Boulesteix et al. (1974) 

3 http://galexgi. gsfc.nasa.gov/Documents/ERO _data_description_2.html 
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images centered on each radio position. Because the far-UV and near-UV images are similar, 
only the near-UV is used in Figures 1-14. All images have an orientation of north up, east 
left. In all figures we have placed a ~9" circle around the center coordinate from Johnson et 
al. (2001) to reflect the positional uncertainties, which are much larger than the positional 
uncertainties of all other images in this study. For all images, the wavelengths are designated 
in the upper right corner. The upper left corners of the near-UV and 8 /im images show 
the object's designation (e.g., 1 for Figure 1). The J and 3.6 fim images give the labeling 
of the objects that we think are the radio-selected sources (discussed below). In the 8 /im 
images we have encircled any extended dust/PAH shells near the radio position. In the 4.5 
/im image, we have superimposed contours of our 6 cm radio maps. 

We identify potential IR counterparts to the radio emission in the Spitzer images. Spitzer 
images are the obvious choice because they are the closest in wavelength to the radio maps. 
Also, if these objects are embedded sources, then the dust emission around these objects 
should be very strong and should obscure emission from shorter wavelengths. Below, we 
discuss the object determination on a case-by-case basis. The celestial positions of these 
objects described below are given in Table 1. 

Object 1: At all wavelengths in Figure 1, excluding UV, the brightest potential source is 
found at the northeast position of the uncertainty circle around the object's position. How- 
ever, the 6 cm emission peaks to the southwest of this object. A source appears coincident 
with this radio emission in the near-UV, 3.6 /im, and 4.5 /im images, and we designate this as 
object 1. Photometry of this radio-coincident source is contaminated by the closer, brighter 
source. The lack of detection at 24 /im of object 1 is probably inconsistent with this object 
being an example of a young, embedded phase of star cluster formation. Object 1 is possibly 
a background object, and therefore, we have eliminated it from further discussion. Johnson 
et al. (2001) state that some of their sample could be background galaxies. 

Object 2: The radio position from Johnson et al. (2001) is located inside a Wolf-Rayet 
nebula, termed BCLMP 638 or IC 132. The Wolf-Rayet star ([MC83] 4 23, a WN star) is 
located in the southwest part of the uncertainty circle in both the 3.6 jum and the B image 
in Johnson et al. (2001). We designate object 2 as the prominent dust/PAH emission part 
of BCLMP 638 located at the southeast portion of the uncertainty circle. The MIPS 24 /im 
emission is bright in this region as well. The other MIPS bandpasses cannot identify whether 
the dust emission arises around the Wolf-Rayet star or object 2 because of lower resolution 
and larger pixel size of the detectors. The 6 cm radio contours suggest that the bright 
radio emission originates around the Wolf-Rayet star and two bright sources, probably stars, 
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9" west of the Wolf-Rayet star. The Ha image indicates that IC 132 is a bipolar nebula. 
Knowing the size and shape of the nebula, we suspect that the bright 6 cm knots of emission 
are found coincidentally on top of the two bright sources in the 4.5 /iin image and that 
the bright 6 cm knots are actually bright radio knots along the edge of the bipolar nebula. 
At JHKs wavelengths, we see the Wolf-Rayet star, but no emission is found at the bright 
dust /PAH emission source to the east of the Wolf-Rayet star. At UV wavelengths, the 
Wolf-Rayet star shows up prominently in the images, but object 2 is not detected at all. 

Object 3: We identify the IRAC object in the middle of the uncertainty circle with the 
radio-selected source. The 6 cm radio contours and JHKs emission are coincident with this 
IRAC source. Some Ha emission is present for this object. Note that the dim IRAC object 
to the southwest of object 3 is just as bright as object 3 at Ha. Johnson et al. (2001) also 
detect this object in their B image, and their radio position is coincident with the H II region 
IC 133 or BCLMP 623. Finally, object 3 is also detected at UV wavelengths, implying little 
extinction for this object. 

Object 4: Like object 1, this object is detected only at IRAC wavelengths and is located at 
the center of the radio uncertainty circle. The 6 cm contours in the 4.5 /xm image show that 
this object is coincident with the radio emission. The lack of detection at 24 /im of object 
4 is probably inconsistent with this object being an example of a young, embedded phase 
of star cluster formation. Object 4 is possibly a background object, and therefore, we have 
eliminated it from further discussion. Johnson et al. (2001) state that some of their sample 
could be background galaxies. 

Object 5: The celestial position of object 5 is coincident with the H II regions BCLMP 
35/36 (Johnson et al. 2001). Two objects are found at the outer edge of the uncertainty 
circle for the radio position. We term the north object "a" and the southern object as 
"b." Looking at the MIPS images, object 5a is the dominant 24 /im source, but object 5b 
is also present but dimmer. The 6 cm contours show that object 5a is the brighter radio 
object but that object 5b has some associated radio emission. The two objects appear to 
be enveloped by a diffuse region of 6 cm emission. Because of this 6 cm connection, we 
have labeled a third object ("5s") as the combination of the two. Object 5s is designated 
in the 8 /xm image. In the JHK S images, we do not detect either of the objects. At Ha, a 
complex ionized gas morphology surrounds these objects. Object 5a is detected at Ha and 
the morphology around this point appears shell-like. Object 5b is the brightest Ha part of a 
shell-like complex south of the uncertainty circle. The UV and B (see Johnson et al. 2001) 
images show a point source in the center of the uncertainty circle that is not associated with 
objects 5a and 5b. This UV-optical point source is close to the center of object 5s. 

Object 6: In the 3.6 /xm image of Figure 6, a star-like source is found at the southeastern 
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portion of the uncertainty circle, and we term this object "6a" . Northeast of the uncertainty 
circle is a source of bright PAH and dust emission, which we term "6b" . In the MIPS 24 
/iin image, these sources are enveloped by diffuse emission and may be interrelated. In the 
other MIPS images, the dominant dust emission appears to be from object 6a. In the 8 fim 
and 24 /xm images, a shell of PAH emission is close to the radio position, which we label as 
object "6s." Faint 6 cm emission surrounds objects 6a and 6b, but no compact emission is 
associated with either IR source. In the JHK S images, we can detect object 6a only. The 
Ha morphology is quite complex with several shells around the radio position. Object 6a 
does not appear to have any Ha emission, while a few points of Ha emission coincident with 
object 6b. Ha emission appears to fill the cavity of object 6s. Object 6a is detected faintly 
in the B image of Johnson et al. (2001), and the celestial position is coincident with the 
H II regions BCLMP 39/40. No UV emission is apparent for any of the objects. 

Object 7: This radio position is coincident with [GDK99] 91, and Gordon et al. (1999) 
claim that this object is a H II region. An IRAC and MIPS point source (object 7a) is found 
almost exactly at the radio position. No 6 cm compact radio emission from our maps is 
associated with this object. The object "7b" is found well north of the radio position and 
is a compact 6 cm source. At 70 /im object 7b dominates the emission. Perhaps object 7b 
is the compact radio source 7 from Johnson et al. (2001). Both objects are present in the 
JHKs images and at Ha. The gas morphologies around those objects appear to be simply 
spherical for both objects. In the GALEX UV images, only object 7b is detected. 

Object 8: In the IRAC images, a bright source appears in the southeastern part of the 
uncertainty circle. No additional sources appear around the radio position in either MIPS or 
IRAC images. This bright IRAC source is also the source of the 6 cm emission. At 8 /im the 
source appears to be associated with a long extended line of PAH/dust emission. However, 
at Ha, this arc of emission disappears, and a simple sphere explains the Ha morphology. 
This source is also detected in the B image in Johnson et al. (2001), and the radio position 
is coincident with the H II region BCLMP 703. Some UV emission is also apparent for this 
object. We did not detect object 8 in the JHKs images. 

Object 9: A bright point source appears in the uncertainty circle at all Spitzer wavelengths 
and in the JHKs images. The 6 cm radio emission is also centered on this Spitzer source. 
A bright Ha source is found on top of the location of object 9. An arc of Ha emission is 
present to the north of object 9. This object is also detected in the B image of Johnson et 
al. (2001), and the radio position is coincident with the H II region BCLMP 79. No UV 
emission is apparent for object 9. 

Object 10: Two point sources appear at opposite edges of the uncertainty circle at 3.6 
fj,m in Figure 10. We have termed the eastern most "10a" and the other "10b." Object 
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10a dominates at 24 /zm, while object 10b dominates the PAH emission at 8 jum. The 6 
cm contours envelope both objects with some radio emission connecting the two points, 
similar to object 5. Thus, we have also labeled a set of flux measurements with "10s" for the 
combination of the two. No sources are detected in the JHKs images. At Ha, both sources 
are on the edge of a shell-like complex of Ha emission. The center of this shell is devoid of 
Ha, but according to the 6 cm contours, some thermal radio emission is apparent. At Ha, 
these objects do not look interconnected as when looking at longer wavelengths. Object 10a 
is detected in the B image in Johnson et al. (2001), and the radio position is coincident with 
the H II region BCLMP 87. Near-UV and far-UV emission emanates from object 10a only 
but not from object 10b. 

Object 11: The only bright star-like source near the radio position is found in the north- 
eastern part of the uncertainty circle in the 3.6 /xm image in Figure 11. We have labeled this 
source as "11a." However, the 6 cm contours of the brightest radio emission are not coin- 
cident with this object but are associated with a ring of dust/PAH emission found around 
the radio position uncertainty circle. This shell is labeled object lis in the 8 jum image. 
The bright dust/PAH sources in this circle are coincident with the 6 cm emission points. 
However, the B image of Johnson et al. (2001) shows a clear source in the middle of the 
uncertainty circle that we see only at UV. This source is not object 11a because we also 
detect object 11a in the UV image. The radio position is coincident with the H II region 
BCLMP 77. 

Object 12: The only source near the radio position is found in the southwestern part of 
the uncertainty circle in the 3.6 fim image of Figure 12. In the MIPS images this source 
appears to be the dominant dust emitter. The 6 cm contours also lie directly on top of 
object 12. Object 12 was not detected in the JHK S images. At Ha object 12 appears to 
have a spherical morphology with no shell-like emission around the object. This object is 
also detected in the B image of Johnson et al. (2001), and the radio position is coincident 
with the H II region BCLMP 714. Object 12 is also detected at UV wavelengths. 

Object 13: In Figure 13 bright 3.6 /im source is found in the middle of the uncertainty 
circle. In the MIPS and IRAC images, no other dominant source exists around the radio 
position. The 6 cm contours are also coincident with the IR source. The Ha emission appears 
spherical for object 13. A crescent-shaped object is located inside the uncertainty circle in 
the northeastern portion of the B image in Johnson et al. (2001), and the radio position is 
coincident with the H II region BCLMP 712. Some UV emission is found to the south of the 
uncertainty circle but is not associated with object 13. 

Object 14: A star-like source appears in the southwestern part of the uncertainty circle 
(see the 3.6 /xm image of Figure 14). We have termed this object "14a." The dust and PAH 
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emission in MIPS and IRAC images appear to be in a cometary V-shape to the south of the 
radio position. Assuming that this V-shaped dust/PAH emission is related, we have labeled 
it "14s" and provided a circle for reference in the 8 /iin image. The 6 cm contours in the 
4.5 /iin image show no significant radio emission around this object. Some faint emission is 
associated with object 14a, but no 6 cm emission appears to be associated with object 14s. 
At Ha the cometary emission of object 14s at 8 /im is only a small portion of the shell of 
Ha emission around object 14a. This shell of emission is oblong with the major axis along 
the northeast. No Ha emission is associated with object 14a. Object 14a is also detected 
in the B image of Johnson et al. (2001), and the radio position is coincident with the H II 
regions BCLMP 749/750. Object 14a is detected at UV wavelengths. 

3.2. Aperture Photometry of the Radio-Selected Objects in M 33 

The aperture used in the near-IR (SQIID and 2MASS), IRAC, Ha, and GALEX UV 
images for the star-like objects (i.e., those that are not shell apertures) had a radius of 
7". The inner and outer radii of the sky annulus for background subtraction were 7" and 
14 " respectively, making the sky area greater than that of the aperture. For object 2, a 
significant bipolar nebula is present in Ha that is not seen with the Spitzer data. For deriving 
a Ha flux, we extended the sky annulus well outside this large nebula to remove only sky 
emission. The shell apertures (i.e., those termed with a "s") along with their sky annulus 
are much larger and set on an individual basis to encompass the brightest emission from 
these shells. Aperture photometry for the "s" designated objects was not done with the 
GALEX UV data, but was done for the Ha images since the shells were detected at Ha. 
The JHK S magnitudes and colors as well as Ha and UV fluxes are presented in Table 1. 
These magnitudes, colors, and fluxes are not corrected for Galactic or internal reddening. 
For those sources not detected in any of the JHKs, Ha, or UV images, we give 3a upper 
limits for the fluxes and lower limits for magnitudes. 

Aperture corrections were mandatory for all IRAC fluxes. We generated our own aper- 
ture corrections for our specific aperture sizes using the theoretical point spread function 
(PSF) images provided by the SSC 5 . Since the aperture and sky annulus were constant for 
all star-like IRAC photometry, the multiplicative values are 1.135, 1.14, 1.176, and 1.285 for 
IRAC 3.6 fim, 4.5 /im, 5.8 /im, and 8 /im respectively. Aperture corrections for the shell 
photometry were also generated with the theoretical PSF images and are ~1. See Table 2 
for the IRAC fluxes. 



5 2004 January In-Flight PSF images downloaded from http://ssc.spitzer.caltech.edu/irac/psf.html 
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For the MIPS images, we used a constant aperture of 48" for all 3 channels. This aperture 
size is the smallest reasonable aperture size for the 160 /im images where the FWHM is ~40". 
Sky subtraction was also performed with the area of the sky annulus greater than or equal 
to the area of the source aperture. The MIPS science images were in the units of DN s _1 . 
The conversion factors from DN s _1 to /iJy arcsec -2 were 1.032, 14,900, and 1,000 for the 24 
jum, 70 /im, and 160 /im images respectively. Aperture corrections were also mandatory for 
all MIPS fluxes. We generated the aperture corrections from theoretical PSF images 6 using 
the same size apertures and sky subtraction areas as those used for the MIPS photometry. 
These aperture corrections were 1.13, 1.61, and 1.89 for the 24 /im, 70 /im, and 160 /im 
images respectively. The MIPS fluxes for the objects are given in Table 2. Because the 
resolution of the MIPS images is lower than those of IRAC and SQIID/2MASS, the MIPS 
fluxes are given in the shell designation for those objects with shells (e.g., lis) or the best 
possible candidate (e.g., object 3). 



3.2.1. Uncertainties 

Uncertainties for all fluxes and magnitudes are given. For the IRAC fluxes, the un- 
certainties account for Poisson noise, variation in the background, correction for correlation 
in the subsample images, and ±6% uncertainty caused by flux calibration. For the MIPS 
fluxes, the uncertainties accounted for are similar to those for IRAC, but the flux calibration 
uncertainty is ±10% for MIPS 24 /im and ±20% for MIPS 70 /im and 160 /im. For the 
JHKs magnitudes, Ha fluxes, and UV fluxes, we account for Poisson noise and variation in 
the background. These uncertainties are presented in Tables 1 and 2. 



3.3. Extinctions, Ages, and Masses Derived From the JHK S Images 

In Figure 15 we plot the colors of the nine radio-selected objects with JHKs detections 
in all three bandpasses. Also plotted are the J-H and H-Ks colors as derived from the 
Starburst 99 v4.0 SEDs (Leitherer et al. 1999). The important parameters for these instan- 
taneous burst models are fixed stellar masses of 10 6 M Q , upper and lower stellar mass limits 
of 100 M Q and 1 M , and a Salpeter initial mass function (IMF). The solid line shows model 
colors for a metallicity of Z = 0.008. The dashed line shows model colors for a metallicity 
of Z = 0.02. The metallicity of each object was estimated using the M 33 oxygen abun- 
dance gradient equation from Garnett et al. (1997). To derive the metallicities, we assumed 



downloaded PSF images from http://ssc.spitzer.caltech.edu/mips/psffits/ on 2004 April 
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that the objects' metallicities relative to solar were proportional to the ratio of their oxygen 
abundances. Since Starburst 99 models have only a few discrete metallicities, we used the 
closest Starburst 99 metallicity, typically Z = 0.008 or Z = 0.02. 

We do not use the colors for object 2 to derive extinctions, ages, and masses because 
the colors found for object 2 are for the Wolf-Rayet star. Thus, the central photoionizing 
source is not a star cluster. We could not determine colors for objects 1, 4, 5, 8, 12, and 13 
because they are not detected in all three near-IR bands. 

3.3.1. Extinctions 

Extinction measurements were estimated for seven objects by extrapolating the mea- 
sured colors back along the reddening vector to the unreddened Starburst 99 model colors. 
We know the cluster's color, the slope of the extinction vector (Schlegel et al. 1998), and the 
theoretical colors. The distance between the cluster color and unattenuated model color was 
converted into a K$ extinction, Ak s (Schlegel et al. 1998). The results of these extinction 
derivations are found in Table 3, and the values range from Ak s of 0.23 mag to 1.31 mag. 
These extinctions, as well as all other derived extinctions, have not had Galactic reddening 
removed. The K s extinction caused by the Galaxy is 0.015 mag approximately (Schlegel 
et al. 1998). Thus, the actual range of reddening caused by the star clusters and M 33 is 
0.21-1.29 mag. 

3.3.2. Ages 

All points on the Starburst 99 model track have an associated age. Moving the observed 
cluster colors back along the reddening vector, the observed colors yield an estimate of the 
cluster age. The la age uncertainties are based from the la color uncertainties. The cluster 
ages and their uncertainties are given in Table 3 and range from 2 Myr to 10 Myr (i.e., young 
clusters) . 

We could not derive an age for cluster 9 because the cluster shows a significant H-K s 
excess indicative of dust emission and/or emission lines. These excesses have been seen in 
star clusters with ages less than 6 Myr (Buckalew et al. 2005). Thus, by comparison to 
previous results, we conclude that the H-K s excess implies an age for cluster 9 of less than 
6 Myr. 
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3.3.3. Masses 

Once the extinctions and ages are known, we can use a similar procedure to the one 
outlined in Buckalew et al. (2005) to derive the masses. Instead of using the blue absolute 
magnitudes, we use unreddened absolute JHK S magnitudes to derive the stellar masses. We 
compared our derived absolute magntiudes with those derived from the 10 6 M Starburst 
99 model at the appropriate age and metallicity. The scaling factor between the model's 
absolute magnitude and the one that we derived gave us the scaling factor to apply to the 
mass. Multiplying the scaling factor to the model mass gives us the derived stellar mass for 
the cluster. Table 3 lists these masses, which have a range of 10 3 -10 35 M . 

3.4. Spectral Energy Distributions 

The UV to mid-IR SEDs are plotted in Figures 16-18. The size of the symbols used 
in the SEDs are equivalent to the flux uncertainties. Upper limits are given as downward- 
pointing arrows. Using these SEDs, we model the free-free and dust emission around these 
objects. 

3.4-1. Dust emission 

Using the 3 MIPS flux results, we derive dust masses and temperatures. We derive 
the dust temperature using simple blackbody fits of the 3 MIPS fluxes. With only 3 data 
points, amplitude and dust temperature are the free parameters. In some cases, the data 
were best fit by two blackbodies with the same amplitude but different temperatures. The 
best fit was found for the set of parameters that minimized the residuals between the model 
and measured fluxes. Two blackbodies were only used when such a fit produced a factor of 
10 decrease in the residuals (objects 2, 3, and 10). The best-fit temperatures are found in 
Table 3 and range between 60 K and 90 K. 

Dust masses were calculated using the procedure outlined in Dale et al. (2001a). If 
only one temperature was used to model the dust emission, the MIPS 70 fim flux was used 
to calculate the dust mass. If two blackbodies provided a better fit, the MIPS 70 /im and 
160 fim fluxes were used to derive the dust masses with the 70 /xm flux used for the hotter of 
the two temperatures. The dust masses are given in Table 3 and range between 1 and 1000 
M . 
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3.5. Extinctions from Ha and UV Fluxes 

Using the Ha fluxes from Table 1 and radio fluxes from Johnson et al. (2001), we 
calculate Ha extinctions using the equations detailed in Caplan & Deharveng (1986). We 
assume a temperature of the H II region of 10,000 K, a N(He + )/N(H + ) of 0.1, densities of 
100 cur 3 , and 100% thermal emission at 6 and 20 cm. We compute extinctions for Ha 
using both the 20 cm and 6 cm fluxes from Johnson et al. (2001). We present the average 
Ha extinctions (A Ha ) in Table 3. The extinctions range from 0.72-3.46 mag. Extrapolating 
these extinctions to Ak s equivalents using the extinction law from Schlegel et al. (1998) 
and removing Galactic reddening, the range of A Ks extinctions is 0.081-0.44 mag. 

We computed UV extinctions by comparing the measured UV emission of the objects to 
GALEX FUV and NUV fluxes calculated from Starburst 99 SEDs. The important param- 
eters (metallicity, stellar mass, and age) for each SED came from the results derived with 
the JHKs colors and magnitudes. The other standard parameters of Starburst 99 were 
assumed (Salpeter IMF, etc.) as they were detailed in §3.3. The GALEX system response 
curves were taken from the GALEX data handbook. The fluxes generated with the filter 
response curves and Starburst 99 SEDs were assumed to be the unreddened FUV and NUV 
fluxes. Ratioing the measured and model fluxes, we calculate extinctions for 0.2271 /im and 
0.1528 /im (A .227i and A i 528 ) and present these extinctions in Table 3. Three objects had 
the necessary UV fluxes and JHK S ages and masses to derive these extinctions. The values 
of A). 2271 are 3.88 mag, 0.686 mag, and 0.15 mag, and the values of A0.1528 are 4.2 mag, 1.9 
mag, and 0.43 mag for objects 3, 11a, and 14a respectively. Converting these to equivalent 
Ak s extinctions using a factor of 15 (based on the Ajj/Ax from Schlegel et al. 1998) and 
removing Galactic reddening, the A Ks range is 0-0.26 mag. 



4. DISCUSSION 

4.1. Morphologies of the Radio-Selected Objects 

These objects in M 33 exhibit a diversity of PAH/dust morphologies in the IRAC 5.8 jum 
and 8 /im images and the MIPS 24 /im image. We combine these IR morphologies with the 
morphologies seen at Ha to determine a multiwavelength classification of the morphologies. 
We find that a simple classification scheme of shell and spherical morphologies applies best 
to these objects. Objects 3, 7, 8, 9, 12, and 13 are classified as spherical, and objects 2, 5, 
6, 10, 11, and 14 are classified as shells. We present an IRAC color-color diagram in Figure 
19 using the 8 /im, 5.8 /im, and 4.5 /im fluxes from Table 2 for all objects. The MIPS fluxes 
from Table 2 are plotted in the color-color diagram in Figure 20. We conclude that no clear 
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separation of the different morphologies can be discerned from either color-color figure. 

4.2. Individual Descriptions of the Objects 

Object 2: A bipolar nebula around a Wolf-Rayet star makes up this thermal radio source. 
We find an Ha extinction (A Ha ) of 0.72 mag. Dust masses and temperatures were 12 M at 
64 K and 330 M at 35 K. The 12 M of dust at 64 K could be ejected stellar material from 
the WN Wolf-Rayet star. The colder dust mass is more likely dust from its natal molecular 
cloud and/or swept-up ISM. 

Object 3: Object 3 (IC 133) is one of the brightest point sources in M 33 at 24 /im. This 
object appears to have a spherical morphology of ionized gas (Ha) and dust (8 /im, 24 /im) 
around the central star cluster. Fitting two blackbodies, we derive dust temperatures of 
87 K and 44 K and dust masses of 31 M Q and 810 M & respectively. The age of this star 
cluster is 7 Myr, and the stellar mass is 10 3 ' 9 M & . The extinction derived from the JHK$ 
color-color diagram is a factor of 1.7 larger than that derived from Ha and the UV. 

Object 5: This multiple-peaked IRAC source is classified as a shell morphology in our 
simplified scheme. From the MIPS data, we derive a dust temperature of 66 K and a mass 
of 12 M Q . The Ha extinction for the shell is 0.51 mag. 

Object 6: Object 6 is interesting for the shell structure seen at 8 /im that lies outside the 
uncertainty circle of the radio position. Object 6a is a stellar source inside the uncertainty 
circle but is found at the northwestern edge of this shell. Object 6b appears as a bright 8 
/im source on the edge of this shell and is also apparent at 3.6 /im. Object 6a is one of the 
youngest star clusters in this sample with an age of 3 Myr, and its derived stellar mass is 
10 3 ' 3 M Q . The dust mass derived for this system is 22 M Q at a temperature of 61 K. The 
derived extinctions from the JHK S and Ha images are similar (A Ks extinctions of 0.23 mag 
and equivalent A Ks of 0.24 mag respectively). 

Object 7: Object 7a is a point source found at the exact location of the radio position. 
However, object 7b is also a bright point source and has bright 6 cm contours on top of its 
position (something not seen with object 7a; see Figure 7). We find that object 7a has a very 
young age of 2 Myr and a stellar mass of 10 3 7 M Q . Object 7b is older at an age of 10 Myr 
and a stellar mass of 10 4 M Q . The dust characteristics, which we derive for the synthesis 
of the two dominant objects, are a dust temperature of 73 K and mass of 7.6 M . Object 
7a is detected at Ha but not UV. Object 7b is detected at both wavelengths. The JHKs 
extinction for object 7a is 1.4 times higher than that derived from Ha (Ax s extinctions of 
0.54 mag and equivalent A Ks of 0.14 mag respectively). The JHK S extinction for object 7b 
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is 3 times higher than the extinction derived from Ha (Ak s of 1.31 mag and equivalent Ak s 
of 0.09 mag respectively). The UV extinction derived for object 7b is 2.5 times smaller than 
that from the JHK S images (A Ks of 0.31 mag and equivalent A Ks of 1.31 mag respectively). 

Object 8: Object 8 is a bright ellipsoid found at the southeastern section of the uncertainty 
circle (see Figure 8). The point source is found along a line of unassociated extended 8 fim 
emission. Because the object was not detected in the JHK$ images, no stellar mass, age, 
or extinctions were derived. We derived an Ah q of 0.36 mag. The dust characteristics are 
typical with a temperature of 66 K and a mass of 16 M Q . 

Object 9: This object is a strong, compact source at all infrared wavelengths and is co- 
incident with the 6 cm radio contours (see Figure 9). The JHK$ color-color diagram in 
Figure 15 shows that this cluster candidate has a significant H-K$ excess. The unusually 
red H— Ks of 1.1 mag may be caused by emission from hot dust or nebular lines in the 
K-band (see Buckalew et al. 2005) and implies an age less than 6 Myr. The J-H of 0.8 mag 
is attributable to extinctions of at least A Ks ~ 0.68 mag. The JHK S extinction is 1.2 times 
higher than the A Ks equivalent extinction derived from Ha. Assuming an age of 6 Myr, we 
find that the stellar mass is 10 4 4 M Q , the highest in the sample. The dust characteristics of 
this object are odd with a dust temperature of 95 K, the hottest in the sample, and mass of 
1.5 M©, the smallest in the sample. 

Object 10: Object 10 is classified as having a shell-like morphology because of its multiple- 
peaked but associated dust/PAH sources, similar to object 5. Because objects 10a, 10b, and 
10s were not detected in the JHK S images, no stellar mass, age, or extinctions were derived. 
The extinctions at Ha were the only derived and had the values of 2.45, 3.46, and 0.96 mag 
for objects 10a, 10b, and 10s respectively. The dust characteristics of this system were best 
fit with two blackbodies with temperatures of 65 K and 40 K and with masses of 55 M Q and 
980 M . 

Object 11: Object 11a is a stellar source inside a shell of dust/PAH emission termed object 
lis. This shell is also the object that is associated with the 6 cm emission (see Figure 11). 
The Ak s of 0.79 mag from the JHKs images is 1.7 times larger than that derived from Ha 
(equivalent A^ s of 0.22 mag). The UV extinctions are 2 times smaller than that derived 
from the JHK S images. The star cluster has an age of 7 Myr and a mass of 10 3 2 M . The 
dust characteristics are a temperature of 57 K and a mass of 64 M Q . The overall morphology 
of this object suggests the presence of a UV-bright star cluster at the center of a wind-blown 
shell, illuminating the surrounding ISM which emits strongly at mid and far-IR wavelengths. 
Based on the age and morphology, this star cluster has similar results to old (an age greater 
than 6 Myr) star clusters from Buckalew et al. (2005). 
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Object 12: Object 12 is a spherical source detected from the UV to mid-IR. We could not 
derive the stellar properties. The extinction from Ha implies a Ak s of 0.21 mag, not a large 
extinction for an embedded source. The dust properties of this source are 15 M at 61 K. 

Object 13: Object 13 is not seen at near-IR nor UV but appears as a point-like source 
at MIPS, IRAC, and Ha. The Ha extinction of object 13 is A Ua = 1.44 mag. The dust 
properties of object 13 are a temperature 58 K and a mass of 23 M Q . 

Object 14: A star-like source (object 14a) appears just outside the southwestern part of the 
uncertainty circle in Figure 14. Object 14a is surrounded by a cometary PAH/dust emission 
and a shell of Ha emission. The age of object 14a is 7.5 Myr, and the stellar mass derived is 
10 3 ' 4 Mq. Such an old star cluster should have a more evolved (i.e., shell-like) morphology 
(see Buckalew et al. 2005). The extinction from JHKs is 1.06 times larger than the Ha 
extinction for object 14a. The UV extinction of object 14a is 1.2 times lower than that from 
the JHKs images. The dust properties of object 14 are a temperature of 56 K and a mass 
of 27 Mr,. 



5. SUMMARY 

These radio-selected objects were thought to be embedded young star clusters, possibly 
larger analogues to ultracompact H II regions. We have observed the positions of these radio- 
selected objects in M 33 from the UV to the mid-IR. We have detected the majority of the 
radio-detected candidates at all of these wavelengths. These 11 radio-detected star clusters 
(excluding objects 1, 2, and 4) are young, when O stars are still present, but appear to have 
ages and extinctions that place them in the detectable limits of optical images. All 11 star 
clusters are known optical H II regions. 

Using the near-IR data, we derived ages, extinctions, and masses. For the limited age 
results, they range from 2 Myr to 10 Myr. The stellar masses derived from the near-IR data 
are typically 10 3 -10 4 M Q . Fluxes from the Spitzer MIPS data were used to generate the dust 
properties of these clusters. We find dust masses ranging from 1-1000 M and temperatures 
between 60 and 90 K. At optical and UV wavelengths, we detected the majority of clusters 
at Ha, 2271 A, and 1528 A. The free-free radii and densities range from 5-20 xlO 17 cm and 
from 1500-9000 e~ cm -3 . Extinctions were derived from the radio, JHKs, Ha, and UV 
fluxes, range from equivalent A Ks values of to 1.29 mag (removing Galactic reddening). 
For an individual object, the maximum difference between various Ak s measurements is a 
factor of 1.6 on average (3 at worst), and considering our assumptions for calculating the 
extinctions and the uncertainties of the data, these factors are acceptable (see §§3.3.1 and 
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3.5 and Table 3). More data and modeling are required to determine what causes these 
variations: geometry or model assumptions. 

The Ha/ Spitzer morphologies show no differences in their IRAC and MIPS colors. Half 
the objects have a young (i.e., spherical) morphology while the others have an older (i.e., 
shell) morphology. Objects 1, 4, 6, 9, and 13 are not detected in the UV, possibly indicating 
larger amounts of extinction than the rest. Objects 1 and 4 are the background objects. 
Objects 9 and 13 are also those with spherical morphologies while object 6 has a shell 
morphology but young (3 Myr) age. 

In general, eleven of the fourteen objects in this radio-selected sample are optically- 
detected star clusters. We must conclude that radio detections of point-like thermal sources 
do not show embedded star clusters, but the brighter free-free emitting H II regions within 
M 33. We propose that infrared images (e.g., MIPS 24 pm) must be used to discern which 
objects in a galaxy are enshrouded in their natal molecular cloud material. These 24 /im 
sources must be checked with optical data (e.g., UBVRI and/or Ha) to determine those that 
are truly embedded. 
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Fig. 1.— Near-UV, Ha, JHK S , IRAC, and MIPS images of object 1. Each image is ~90" on 
a side with an orientation of North up, and east left. The radio position from Johnson et 
al. (2001) is designated by the circle in the center of the images. The radius of the circle 
reflects roughly the uncertainty in their position. In the 4.5 /xm image, we have overlaid the 
6 cm map from Duric et al. (1993). The candidates are labeled in both the J and 3.6 jum 
images. We located object 1 as the point source on the upper left of the circle. However, 
the 6 cm radio source is found to that candidate's lower right. A very faint source does exist 
closer to this position but is detected only at near-UV, 3.6 /iin, and 4.5 /im. We think that 
object 1 is a background object, possibly an AGN or QSO. 
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Fig. 2. — Near-UV, Ha, JHKs, IRAC, and MIPS images of object 2. Orientation and 
notation is the same as Figure 1. Object 2 is located in the lower left of the radio position 
uncertainty circle. The Wolf-Rayet star is found at the western edge of the uncertainty circle. 
The 6 cm radio contours in the 4.5 /xm image arise on the western side of the bipolar (see 
Ha) nebula. No radio emission appears to be associated with the bright IRAC object. 
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Fig. 3. — Near-UV, Ha, JHK$, IRAC, and MIPS images of object 3. Orientation and nota- 
tion is the same as Figure 1. The object is found in the uncertainty circle at all wavelengths 
and is associated with the 6 cm point source. 
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Fig. 4. — Near-UV, Ha, JHK S , IRAC, and MIPS images of object 4. Orientation and 
notation is the same as Figure 1. The object is found nearly in the center of the radio position 
uncertainty circle in the IRAC images only. The 6 cm emission contours are centered on this 
object in the 4.5 /im image. We think that it is a background object, possibly an AGN or 
QSO. 
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Fig. 5. — Near-UV, Ha, JHKs, IRAC, and MIPS images of object 5. Orientation and 
notation is the same as Figure 1. Two objects, 5a and 5b, appear on opposite sides of the 
uncertainty circle. Because the free-free emission contours connect both objects, we think 
that both are associated, and the combination of the two gives rise to the radio emission. 
This combination of the two objects is termed object 5s and is labeled in the 8 /iin image. 
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Fig. 6. — Near-UV, Ha, JHKs, IRAC, and MIPS images of object 6. Orientation and 
notation is the same as Figure 1. A shell morphology is seen at 8 /im and 24 /iin and 
designated 6s. Object 6b is the brightest PAH/dust emitter, but a stellar source (object 6a) 
is found within the uncertainty circle. No 6 cm emission appears to be centrally concentrated 
in the middle of the image. 
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Fig. 7.— Near-UV, Ha, JHK S , IRAC, and MIPS images of object 7. Orientation and 
notation is the same as Figure 1. Object 7a is the actual star cluster, but no 6 cm emission 
appears to be associated with this object. Object 7b is the only nearby 6 cm emitter and 
will contaminate MIPS 70 /im and 160 /xm fluxes of object 7a. 
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Fig. 8. — Near-UV, Ha, JHK S , IRAC, and MIPS images of object 8. Orientation and 
notation is the same as Figure 1. The only cluster candidate is found within the radio 
position uncertainty circle, and the 6 cm radio contours are exactly located upon this object. 
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Fig. 9. — Near-UV, Ha, JHK S , IRAC, and MIPS images of object 9. Orientation and 
notation is the same as Figure 1. The only cluster candidate is found within the radio 
position uncertainty circle, and the 6 cm radio contours are exactly located on top of the 
IRAC object. 
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Fig. 10.— Near-UV, Ha, JHK S , IRAC, and MIPS images of object 10. Orientation and 
notation is the same as Figure 1. The large circle in the 8 //m image designates object 10s. 
Since the two possible objects appear to be connected via their 6 cm emission, we think that 
the combination of the two (termed 10s) is the actual object. 
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Fig. 11. — Near-UV, Ha, JHK S , IRAC, and MIPS images of object 11. Orientation and 
notation is the same as Figure 1. The large circle in the 8 /xm image designates object lis. 
Object 11a is the stellar source, but the PAH/dust emission for this cluster is emanating 
from a shell around object 11a. The 6 cm contours superimposed on the 4.5 /im image 
suggest that the free- free emission emanates from the shell of dust/PAH emission around 
object 11a. 
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Fig. 12.— Near-UV, Ha, JHK S , IRAC, and MIPS images of object 12. Orientation and 
notation is the same as Figure 1. Object 12 is found at the southwestern part of the uncer- 
tainty circle, is the source of the brightest 6 cm emission, and appears to have a spherical 
morphology. 
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Fig. 13.— Near-UV, Ha, JHK S , IRAC, and MIPS images of object 13. Orientation and 
notation is the same as Figure I. This object with a spherical morphology is found at the 
center of the radio position uncertainty circle and is the source of the brightest 6 cm emission. 
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Fig. 14. — Near-UV, Ha, JHK S , IRAC, and MIPS images of object 14. Orientation and 
notation is the same as Figure 1. The large circle in the 8 /xm image designates object 14s. 
While 14a is the stellar component, the PAH and dust emission appears to be coming from 
the cometary-like emission around the radio position (termed 14s). Some 6 cm emission 
appears to be associated with the stellar source, but no emission is associated with the 
cometary PAH/dust emission. 



-35 - 



Fig. 15. — JHK S color-color diagram of eight objects. The legend at the upper left indicates 
which symbol goes with which object. Uncertainties are la in length for all colors. The 
connected crosses show the evolutionary model tracks for instantaneous star forming bursts 
from Starburst 99 for metallicity Z = 0.008 (solid) and for metallicity Z = 0.02 (dotted). 
Crosses along the model tracks indicate ages in step of 1 Myr. Specific ages are labeled along 
the tracks (0, 6, 10, and 20 Myr). The arrow to the right shows the reddening vector for 

Ay = 5. 
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Fig. 16. — Spectral energy distributions of objects 2, 3, 5, and 6 from the ultraviolet to 
the mid-infrared. Uncertainties for these fluxes are the size of or smaller than the symbols 
plotted. Downward-pointing arrows are upper limits. If multiple objects are designated for 
one position, a legend is provided to indicate which symbol goes with which object. 
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Fig. 17. — Spectral energy distributions of objects 7, 8, 9, and 10 from the ultraviolet to the 
mid-infrared. Notation is the same as that for Figure 16. 
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Fig. 18. — Spectral energy distribution of objects 11, 12, 13, and 14 from the ultraviolet to 
the mid-infrared. Notation is the same as that for Figure 16. 
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Fig. 19. — IRAC color-color diagram of the M 33 radio-selected objects. The different 
symbols represent the morphologies discussed in §4.1. Note that no clear separation occurs 
in the IRAC colors for the different morphologies. 
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Fig. 20. — MIPS color-color diagram of twelve of the M 33 objects. The point-by-point 
labeling is the same as that of Figure 19. Note that the morphologies discussed in §4.1 are 
not separable by their MIPS colors. 



Table 1. Radio- Selected Object Positions, JHK S Magnitudes and Colors, and Ha and UV Fluxes 



Object Name 


R. A. (J2000) 


Dccl. (J2000) 


J (mag) 


H (mag) 


K s (mag) 


J-H (mag) 


H-K S (mag) 


F„(Ha) [lO" 5 Jy] 


F„(0.1528 /jm) [10~ 5 Jy] 


F„ (0.2271 |im) [10~ 5 Jy] 


1 


1 33 02.6 


30 46 46.1 














0.85±1.73 


0.772±2.9 


2 


1 33 16.1 


30 56 42.4 


16.33±0.14 


16.3±0.19 


16.19±0.24 


0.028±0.235 


0.114±0.308 


8.27±0.03 


78.4±16.7 


67±27 


3 


1 33 16.6 


30 52 49.9 


15.65±0.1 


14.81±0.08 


14.278±0.08 


0.847±0.124 


0.528±0.112 


17.8±0.01 


23.3±9.1 


21.8±15.4 


4 


1 33 37.5 


30 47 18.5 


>17.5 


>17.1 


>17.6 






<0.8 


<0.08 


<0.05 


5a 


1 33 39.1 


30 38 08 


>17.3 


> 17.6 


> 17 






1.92±0.01 






5b 


1 33 39.3 


30 38 02.2 


>16.5 


>16 


>15.5 






4.47±0.01 






5s 


1 33 39.2 


30 38 03.4 


>16.5 


>16 


>15.5 






12.7±0.02 


62.8±14.9 


68±27 


6a 


1 33 43.6 


30 39 05 


16.43±0.06 


16.12±0.18 


15.61±0.12 


0.31±0.19 


0.51±0.21 


3.68±0.02 


<0.4 


<0.3 


6b 


1 33 44 


30 39 13.5 


>18.6 


>16.9 


>15.3 






<0.5 


<0.4 


<0.3 


6s 


1 33 44.1 


30 39 03.5 


>14.4 


>16.9 


>16 






14.2±0.05 






7a 


1 33 48.1 


30 39 18.2 


16.49±0.12 


15.86±0.2 


15.04±0.1 


0.63±0.23 


0.82±0.22 


2.7±0.008 


<0.3 


<0.2 


7b 


1 33 48.3 


30 39 35.5 


18.54±0.96 


16.5±0.34 


15.5±0.15 


2.04±1.02 


1±0.37 


3.85±0.01 


6.6±4.9 


4.94±7.34 


8 


1 33 59.8 


30 32 44.2 


>17.4 


>16.9 


>17.9 






5.15±0.007 


6.52±4.8 


7.81±9.23 


9 


1 34 00.2 


30 40 47.5 


16.16±0.07 


15.4±0.08 


14.32±0.03 


0.76±0.1 


1.08±0.08 


4.67±0.01 


<0.2 


<0.2 


10a 


1 34 02.6 


30 38 38.7 


>16 


>13.8 


>14 






6.9±0.04 


23.7±9.2 


25.7±16.7 


10b 


1 34 02 


30 38 39 


>16 


>13.8 


>14 






2.73±0.02 


<0.3 


<0.3 


Ids 


1 34 02.2 


30 38 40.8 


>16 


>13.8 


>14 






27.4±0.05 






11a 


1 34 06.6 


30 41 47.3 


15.73±0.05 


15.98±0.27 


15.38±0.07 


-0.25±0.28 


0.6±0.28 


4.98±0.02 


102±19 


79.9±29.5 


lis 


1 34 06.6 


30 41 47.3 


>15.3 


>15.6 


>14.1 






27.9±0.02 






12 


1 34 13.5 


30 34 49.3 


>17.3 


>15.8 


>15.3 






5.25±0.01 


7.54±5.17 


8.37±9.55 


13 


1 34 17.5 


30 33 45.8 


>17.4 


>16.4 


>16 






4.6±0.01 


<0.6 


<0.4 


14a 


1 34 38.6 


30 43 56.8 


15.49±0.02 


15.01±0.03 


14.82±0.02 


0.479±0.032 


0.187±0.033 


3.22±0.01 


72.1±16 


h£fi±28 


14s 


1 34 39 


30 43 55.1 


>15 


>14.6 


>14.5 






8.21±0.01 







Note. — Units of right ascension are hours, minutes, and seconds, and the units of declination arc degrees, arcminutes, and arcseconds. All magnitudes and fluxes have not been dcrcddcncd. Numbers after the ± 
arc lex uncertainties. 



-42- 



Table 2. IRAC and MIPS Fluxes of M 33 Radio- Selected Objects 



Object Name 


IRAC Ch. 1 
3.6 ^m (mjy) 


IRAC Ch. 2 
4.5 (im (mjy) 


IRAC Ch. 3 
5.8 fim (mjy) 


IRAC Ch. 4 
8 (im (mjy) 


MIPS Ch. 1 
24 (im (mjy) 


MIPS Ch. 2 
70 /im (mjy) 


MIPS Ch. 3 
160 (im (mjy) 


2 


0.383±0.032 


0.374±0.031 


1.16±0.1 


3.2±0.3 


57±1 


1190±270 


940±200 


3 


6.22±0.51 


11.6±1 


32.8±2.7 


92.7±7.6 


1700±200 


7340±1810 


4150±880 


4 


0.182±0.016 


0.228±0.02 


0.221±0.02 


0.259±0.023 


<0.08 


<5 


<200 


5a 


0.662±0.055 


0.557±0.047 


4.03±0.34 


12.1±1 








5b 


0.996±0.083 


0.653±0.055 


4.13±0.34 


11.6±1 








5s 


2.97±0.24 


2.01±0.17 


13.7±1.4 


39.7±3.2 


142±21 


1310±320 


1360±340 


6a 


0.155±0.016 


0.152±0.014 


0.293±0.037 


0.981±0.11 








6b 


0.529±0.045 


0.417±0.035 


2.11±0.18 


6.17±0.51 








6s 


3.54±0.3 


2.66±0.23 


17.7±1.4 


51.4±4.2 


101±10 


1810±450 


2080±480 


7a 


0.79±0.066 


0.508±0.04 


2.59±0.21 


8.11±0.68 








7b 


1.53±0.13 


1.24±0.1 


8.17±0.68 


24.1±2 


194±43 


1120±280 


1210±320 


8 


0.749±0.062 


0.585±0.049 


3.44±0.28 


9.75±0.8 


42.1±7.9 


1730±370 


2550±129 


9 


5.85±0.48 


6.12±0.51 


31.7±2.7 


99.4±0.82 


829±182 


426±95 


3230±640 


10a 


0.825±0.069 


0.691±0.058 


4.59±0.38 


12.6±1 








10b 


1.58±0.13 


1.29±0.11 


8.24±0.68 


24.6±2 








10s 


6.08±0.51 


4.97±0.41 


34±3 


93.2±7.6 


254±29 


5560±1050 


3960±850 


11a 


0.32±0.027 


0.339±0.028 


<0.04 


<0.1 








lis 


2.74±0.23 


2.12±0.17 


8.6±0.7 


24.4±2 


75.1±7.9 


4120±1000 


1790±410 


12 


1.04±0.08 


0.852±0.07 


4.07±0.34 


11±0.8 


108±11 


1260±260 


3470±750 


13 


1±0.08 


0.745±0.062 


4.72±0.39 


12.9±1.1 


93.9±17 


1580±390 


3600±780 


14a 


0.529±0.044 


0.318±0.027 


0.908±0.076 


2.3±0.2 








14s 


1.02±0.08 


0.654±0.055 


2.94±0.24 


7.42±0.61 


41.1±5.5 


1620±260 


1280±310 



Table 3. Extinctions derived from JHK S , Ha, UV, and Radio Data along with Derived Stellar Ages, Stellar Masses, 

Dust Temperatures, and Dust Masses 



Object Name 


A Ks (mag) 


Au a (mag) 


Ao.2271 (mag) 


Ao.1528 (mag) 


Stellar Age (Myr) 


Log(Stellar Mass) [Log(M )] 


Dust Temperature (K) 


Dust Mass (M ) 


12 




0.72 






7 +l 

'-2 




64(35) 


12(330) 


3 


0.84 


2.01 


3.88 


4.2 


3.5 


87(44) 


31(810) 


45a 




2.57 














5b 




1.65 














5s 




0.51 










66 


12 


6a 


0.23 


1.8 






o + l 
6 -3 








6b 


















6s 




0.33 










61 


22 


7a 


0.54 


1.06 






2±2 


3.15 






7b 


1.31 


0.68 


4.18 


4.6 


1n +l000 

±u -l() 


3.14 


73 


7.6 


8 




0.36 








66 


16 


9 


0.68 


3.4 






<6 


3.4 


92 


1.5 


10a 




2.45 












... 


10b 




3.46 












^ ... 


10s 




0.96 






7 +2 
'-3 




65(40) 


^55(980) 


11a 


0.79 


1.65 


0.686 


1.9 


3.2 






lis 















57 


64 


12 




1.6 










61 


15 


13 




1.44 










58 


23 


14a 


0.3 


1.76 


0.15 


0.43 


- c+1000 
'•°-0.5 


3.4 






14s 




0.75 










56 


27 



Note. — Numbers after the ± are Icr uncertainties. Extinctions recorded in this Tabic do not have the Galactic reddening component removed (an Ak s of 0.015 mag). Ak s is derived 
from the distance between the Starburst 99 model JHKg colors and the reddened JHKg colors. Ajj a is derived from the radio and Ha fluxes used in the equation from Caplan & 
Deharvcng (1986). Ao.i528/tm and Ao.227l/im were derived by comparing Starburst 99 fluxes at these wavelengths with those measured (see Table 1 for UV fluxes). Dust temperature and 
dust masses found in the parentheses are the second results from the two blackbody fits. 



This figure "fla.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0509886vl 



This figure "f lb.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0509886vl 



This figure "flc.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0509886vl 



This figure "f2a.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0509886vl 



This figure "f2b.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0509886vl 



This figure "f2c.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0509886vl 



This figure "f3a.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0509886vl 



This figure "f3b.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0509886vl 



This figure "f3c.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0509886vl 



This figure "f4a.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0509886vl 



This figure "f4b.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0509886vl 



This figure "f4c.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0509886vl 



This figure "f5a.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0509886vl 



This figure "f5b.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0509886vl 



This figure "f5c.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0509886vl 



This figure "f6a.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0509886vl 



This figure "f6b.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0509886vl 



This figure "f6c.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0509886vl 



This figure "f7a.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0509886vl 



This figure "f7b.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0509886vl 



This figure "f7c.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0509886vl 



This figure "f8a.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0509886vl 



This figure "f8b.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0509886vl 



This figure "f8c.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0509886vl 



This figure "f9a.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0509886vl 



This figure "f9b.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0509886vl 



This figure "f9c.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0509886vl 



This figure "flOa.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0509886vl 



This figure "flOb.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0509886vl 



This figure "flOc.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0509886vl 



This figure "flla.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0509886vl 



This figure "fllb.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0509886vl 



This figure "fllc.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0509886vl 



This figure "fl2a.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0509886vl 



This figure "fl2b.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0509886vl 



This figure "fl2c.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0509886vl 



This figure "fl3a.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0509886vl 



This figure "fl3b.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0509886vl 



This figure "fl3c.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0509886vl 



This figure "fl4a.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0509886vl 



This figure "fl4b.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0509886vl 



This figure "fl4c.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0509886vl 



This figure "fl5.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0509886vl 



This figure "fl6.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0509886vl 



This figure "fl7.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0509886vl 



This figure "fl8.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0509886vl 



This figure "fl9.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0509886vl 



This figure "f20.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0509886vl 



